Diamagnetic response of Aharonov-Bohm rings: impurity backward scatterings 
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We report a theoretical calculation on the persistent currents of disordered normal-metal rings. 
It is shown that the diamagnetic responses of the rings in the vicinity of the zero magnetic field 
are attributed to multiple backward scatterings off the impurities. We observe the transition from 
the paramagnetic response to the diamagnetic one as the strength of disorder grows using both the 
analytic calculation and the numerical exact diagonalization. 
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I. INTRODUCTION 

The phase coherence of electrons in mesoscopic rings 
gives rise to many interesting quantum phenomena. 
Among them the persistent current in clean rings in- 
duced by the Aharonov-Bohm (AB) flux has been ex- 
tensively studied and is now well understood both the- 
oretically and experimentally. However, recently, new 
puzzling behaviors have been observed for the aver- 
aged persistent currents in diffusive rings. Experiments 
clearly indicate that the averaged persistent currents 
show diamagnetism'sr^.'S-: in the vicinity of the zero mag- 
netic field in contrary to the existing theory and the am- 
plitude of the persistent currenli^ is much larger than the 
theoretically predicted values (S^i 

There have been several attempts to explain these puz- 
zling behaviors. The repulsive electron-electron interac- 
tions may enhance the amplitude of the persistent cur- 
rent, but produce only paramagnetic responses for the 
ensemble averaged currents The nonequilibrium ac 
noise in a mesoscopic ring can cause both the decoher- 
ence of electron wave function and the diamagnetic dc 
current in the loopji2i Including the spin-orbit scattering 
in a diffusive ring with ac noise, it was claimed that the 
sign of persistent current changes from diamagnetic to 
paramagnetic as the strength of the spin-orbit scattering 
increasesiii However experiments show that the persis- 
tent currents of Au rings^ and Ag rings^ with strong spin- 
orbit interactions as well as of GaAs rings^ with a weak 
one all exhibit the diamagnetic responses. Even though 
the experiments were performed at much higher temper- 
ature than the superconducting transition temperatures, 
it was claimed that the superconducting fluctuations due 
to the phonon mediated attractive interaction might lead 
to a diamagnetic magnetization.^^ Furthermore an at- 
tempt which included the large contribution of far levels 
from the Fermi level-^ claimed that superconducting fluc- 
tuations would give a much larger persistent current. 

The magnetic response of an AB ring depends on the 
parity of the number of electrons n in the ring. For 
n = AN the states of the topmost electrons, the clock- 
wise moving and the counterclockwise moving states, are 
degenerate when AB flux $cxt is zero. The paramagnetic 



current emerges when the degeneracy becomes lifted due 
to a finite $ext and, thus, the circulating direction of the 
topmost electrons is determined. The ensemble average 
of the currents of rings with n — AN, n = AN ± 2 and 
n = AN ± 1 show paramagnetic responses with the peri- 
odicity $o/2 in clean limit or in weakly disordered case, 
where $0 = h/e is the unit flux quantum. However, in a 
disordered ring, the backward scattering of electrons off 
the impurities may induce transitions between the two 
states and thus severely reduce the paramagnetic persis- 
tent current for n = AN. 

In the disordered systems the backward scattering pro- 
cess occurs through multiple scattering processi^ which 
leads to rich interesting physical phenomena such as weak 
localization^^ and universal conductance fluctuation.— 
Experiments for diffusive normal-metal ringsii as well 
as wiresi^ showed that the resistance oscillation along 
with the external magnetic fleld can be explained in 
terms of the weak localization theoryiiS Therefore we can 
think that there is much higher backscattering proba- 
bility in those normal-metal rings. A numerical study 
showed that the amplitude of backscattering in a disor- 
dered multi-channel wire increases along with the length 
of the wire and the reflection coefficient grows up to be- 
come comparable to the transmission coefficient.^ 

In present study we show that the backscattering pro- 
cess in the diffusive ring induce the transition between 
the topmost levels and, as a result, the ensemble averaged 
persistent current changes from paramagnetic to diamag- 
netic as the impurity strength grows. 



II. DIAMAGNETISM DUE TO 
BACKSCATTERING PROCESS 

The Hamiltonian of an AB ring with potential impurity 
scattering can be given by 
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where km — 2'K(rn + f)/ L and p — 2iim' / L with / = 
*i'ext/*I*Oi the circumference of the ring L and integers m 
and m'. Here Vp is the Fourier component of the impurity 
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potential v(x). The scattering term of the Hamiltonian, such as 

p ^p^\ +p'^km J be divided into the forward scat- 
tering and the backscattering process depending on the 
value of the momentum transfer p. For a diffusive AB 
ring the forward scattering may reduce the amplitude of 
persistent current but cannot transfer electrons to oppo- 
sitely moving states. Thus, although we will consider all 
the scattering processes in the following exact diagonal- 
ization calculation, we first consider only the backward 
scattering process for transparency of argument. 

The backscattering process in an AB ring can be de- 
scribed by cLj,^^ Cfc^, where q = 27r(n + 2f)/L with 
integer n. When there is no threading AB flux and thus 
q — 21m/ L, it is sufficient to consider the most dominant 
process such as cLj, Ck,^ for g = with n = at ground 
state. However, since the value of q cannot be zero for 
finite AB flux due to the symmetry breaking of the mo- 
mentum levels, we need to consider the backscattering 
process with finite q. 

For a multichannel wire the backscattering processes 
take place through multiple scatterings between the 
states in different channels, which accompany momentum 
transfer between the initial and the final state. When an 
electron becomes backscattered eventually, the effective 
backscattering process can be given by w{q)i^_^ +gCfc^^. 
If there is no AB flux threading the loop, the effective 
backscattering can be described as a;(g)cl^ Cfc^ approx- 
imately with ijj{q) = ujQdqfi. However, for finite $oxtj ^{<l) 
will have some distribution which is centered about g « 0. 
For example, the backscattering processes between edge 
states in quantum Hall system through multiple scatter- 
ings off impurities show a Lorentzian distribution^ and 
light scatterings from a one-dimensional rough random 
metal a Gaussian-type distribution for the direction of 
the reflected waves.-- 

Here we take a Gaussian- type distribution for 




FIG. 1: Energy levels of the AB ring. When n = 4,N , (a) 
shows the levels for / ~ and (b) for / ~ 0.5. (c) and (d) 
correspond to the case n = AN — 2. Here filled (open) circles 
denote the occupied (unoccupied) levels. 



-aq a 



(2) 



where a is the lattice constant of the loop and we set 
a = 0.1 here and after. Different values of a or different 
distribution functions result in only quantitatively differ- 
ent behaviors. Then the effective single channel Hamilto- 
nian where the backscattering term comes from multiple 
scatterings between different channels can be written as 
follows, 



km ,(3 
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where = £0(111 + /)^ and Eq = 27r^?i^/r7ieL^ with the 
mass of an electron mg. 

Since the scatterings off the impurities occur mainly 
for the electrons near the Fermi level, it is sufficient to 
consider the scattering processes between several levels 
near the Fermi level as shown in Fig. ^ The contribution 
of other electrons to the persistent current are obtained 
assuming that they do not participate in the scattering 
processes. 

When an electron is backscattered from the state, \km), 
to the state, \km'), the scattering matrix element be- 
comes Wmm' = ^{<l) with km' = —km + 9- The Green's 
function Gmm{^) can be obtained from the equation 
(e ± is — H)G = I which can be written as follows. 
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In Fig. n^a), we consider scatterings between four levels, 
TOi = —N + 1, 7712 = "-^i ™3 — and 7774 = iV — 1, 
with a small external flux, / « 0. Including the electrons 
at the zero momentum state, we can see that the total 
number of electrons 77 in Fig. ^a) is 47V. The scattering 
processes are expressed as arrows in Fig. ^a). 

The Green's functions can be obtained from Eq. Q 
as follows, 
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The Green's functions Gm-smsi^) and Gmimi(e) can be 
obtained from Eqs. Q and © by exchanging the mo- 
mentum indices; m2 ^ and ^ mi. 

All these four Green's functions have the same denom- 
inator 2? represented as 



V — T T 



-(e ±is- e™J(e ±is- tmi)W, 
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from which we can obtain four poles e'^ such that < 
< < e*. The persistent currents are calculated 
using the residues of the Green's functions obtained by 
z^. = \im^^^k{e — e^)Gmimi{^)- Since we consider three 
scattering particles in four states of Fig. Cta), the prob- 
ability that a certain level, m^, is occupied at the ground 
state is represented as the sum of the residues of the three 
lowest energy states, Prm — 



fc=l ^m, 



with Up ~ i. 



Since electrons which are not shown in the figure are 
assumed not to participate in the scattering process, the 
contributions of these electrons can be obtained simply 
summing the momentum of each electron^ 

= Et^"^ + (11) 



with Jo = evp/L and vp = 2-kN/L. Therefore, the to- 
tal persistent current is written as the sum of the two 
contributions as follows, 



' 4 

.4=1 



(12) 



where the first part represents the persistent current of 
the three scattered electrons in Fig. ^a) and the second 
one that of the unscattered electrons from m/ = — A'' -I- 2 
to rrir = N — 2. The prefactor 2 comes from the spin 
degeneracy. 

When the AB flux / increases to / w 0.5, the energy 
levels become shifted as shown in Fig. ^b). In this case, 
since the energy levels of the two states, | A^) and | — A^ — 
1), are much closer than before, the scattering process 
between these two states becomes important unlike to the 
case of Fig. ^a). For the persistent current of Fig. Q^b), 
we use the definition of parameters; mi — —N + 1, m2 — 
-N, ma = TV, m4 = A - 1, m; = ~N + 1, rrir ^ N - 2 
and Up = 2. Here and after we set A^ = 100. 

Since we consider a quadratic dispersion, e.,„ — E^{m + 
fY, the level spacing at the Fermi level is of the order 
of A^£^o- Here we define the level spacing A = NE^. In 
Fig. 121 we show the persistent current /47V for n — AN 
with wo/A — ujq/NEq — 0.7 and a — 0.1, where the left 
(right) curve corresponds to Fig. ^a) (Fig. db)). In the 
intermediate regime of the external flux /, we extend the 
results from the both limits to obtain the approximate 
results as shown in Fig. [3 
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FIG. 2: Persistent currents for a = 0.1 and ljo/A = 0.7. Left 
(right) curves in /4jv and I4N-2 correspond to Fig. Qa) and 
Fig. 0^c) (Fig. H^b) and Fig. H^d)), respectively. We show 
4:Iav for the ensemble averaged persistent current instead of 
lav for clarity. 



When the total number of electrons is n = AN — 2, the 
persistent current I4N-2 can also be obtained similarly 
using the scattering processes shown in Fig. and 
(d), where we use the parameters mi = — A + l,m,2 ~ 
-~N,m3 = A", TO4 = N-1, mi = -A^+2,m,. = A^-2 and 
Up — 2 for Fig. n^c) and mi = —A + 1, m2 = —A, ma — 
A^ - 1, m4 = A - 2, m,; = -A + 2,m,. = A^ - 3 and 
Up = 3 for Fig. IHd). 

The dotted line in Fig. |21 shows the persistent current 
/4Ar_i for n = AN — 1 which is obtained by the rela- 
tion liN-i — 0.5(/4Ar -f- /4Ar_2)- This relation can be 
easily understood for clean loops without impurities. If 
we consider one additional electron at m = — A state in 
Fig. n^c), say spin- up electron, it becomes the configura- 
tion for n = AN — 1 with the number of spin- up electrons 
n„ = 2 A and spin-down electrons Ud = 2N —1. The con- 
figuration for the spin-up electrons is the same as that of 
Fig. n^a) with half the number of electrons and, for the 
spin-down electrons, it corresponds to Fig. ^c). Since 
we do not have spin flip processes, the total persistent 
current is equal to the sum of half the persistent currents 
of n = 4 A and n — AN — 2. In the present case of loops 
with impurities, if the impurity scattering is not spin de- 
pendent so that each spin degree of freedom is decoupled 
from the other, we have the same relation. 

The current /4Ar+i for n = AN -I- 1 is the same as that 
of /4Ar_i with small difference of 0{1/N). Hence, the 
ensemble averaged persistent current lav can be written 
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The ensemble averaged persistent current in clean rings 
shows the paramagnetic responses for / « 0, since the 
paramagnetic currents of I^pf dominate over the diamag- 
netic current Ian -2- Even for an AB loop with impu- 
rities, the qualitative picture does not change although 
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the scattering is expected to reduce the magnitude of 
the paramagnetic current, if one considers only the for- 
ward scatterings off the impurities. However, when we 
consider contributions from the backward scattering pro- 
cesses, the situation changes drastically as indicated in 
Fig. El 

It is well known that the paramagnetic response for 
n = AN comes from the contribution of the electrons with 
TO2 = - A^ in Fig. [H^a). At / = the levels at m2 = -N 
and ms — N in Fig. da) are degenerate. Finite external 
flux lifts the degeneracy and the electrons at ni2 = —N 
contribute a large magnitude of persistent current in ad- 
dition to the currents by the other electrons in the lower 
levels, thus resulting in the paramagnetic ensemble av- 
eraged current. However, when the backward impurity 
scattering is sufficiently strong, the impurity scattering 
causes mixing of the two levels. As a result, the popula- 
tion probability of the level m2 = —N decreases, whereas 
that of the level ms = N increases. This reduces the 
paramagnetic I^n and results in the diamagnetic aver- 
aged current lav for small / in Fig. |21 



III. NUMERICAL CALCULATION 

Since we considered impurity scatterings of only four 
levels in calculating the Green's functions, the persistent 
current curves in Fig. El are only approximate. In or- 
der to obtain more accurate persistent current curves we 
perform a numerical calculation using the exact diago- 
nalization method, where we include the forward scat- 
terings, X]fe„ q '"^(9)ci^+5Cfe^, as well as the backscatter- 
ings. We consider a many-body basis with twelve levels 
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FIG. 3: Momentum distribution of the AB loop for n = 
with 12 levels of which 7 levels are filled when / = 0.1. The 
horizontal axis denote the absolute value of the wave vector 
in the loop. The filled (open) marks show the occupation 
probability of the levels with negative (positive) values of wave 
vectors. Inset shows the case for n — AN — 2 with 12 levels of 
which 6 levels are filled. 
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FIG. 4: Persistent currents for (a) ljo/A = 1.2, (b) t^o/A = 
1.6 and (c) ljo/A = 2.0 obtained by the exact diagonalization 
method with 12 levels. Ian-\ (the dotted line) is obtained 
using the relation I/^n-\ = 0.5(/4jv 4- I/^n-2)- 



near the Fermi level, where seven (six) levels are filled 
for n — AN [n = AN — 2). If there is no impurity, the 
ground state is |111111100000) (|111111000000)) in the 
occupation- number space for n — AN {n — AN — 2). Here 
we does not consider spin degree of freedom of electrons 
and the amplitude of the persistent current will be mul- 
tiplied by the spin degeneracy at the end of calculation. 

We calculate the ground state by diagonalizing the 
Hamiltonian matrix with 792 (924) basis for several val- 
ues of uJo for n = AN [n = 4A^ — 2) and show the calcu- 
lated momentum distribution for n = AN as a function 
of the absolute value of momentum level when / = 0.1 
in Fig. The filled marks denote the levels with neg- 
ative momentum (left branch) and the open marks with 
positive momentum (right branch). For weak impurity 
strength, we note that the lowest seven levels are almost 
occupied. The seven levels are composed of four levels in 
the left branch and three levels in the right branch. In 
this case, it is easily seen that the contribution from the 
topmost level with ka/ir — —0.978 in left branch becomes 
dominant. 

In Fig. Ela), we show a paramagnetic persistent cur- 
rent for weak impurity scattering strength cjo/A = 1.2, 
where the paramagnetic persistent current Inq dominates 
over the diamagnetic one Here, the contributions 

of the electrons which do not participate in the scat- 
tering processes to the persistent currents are included 
as in the second term of Eq. H12|) . However, as the 
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strength of the impurity increases, the occupation prob- 
abihty of the topmost level in the left branch decreases 
while that of the topmost level with ka/ir — 0.982 in the 
right branch increases as can be seen in Fig. |21 The 
persistent currents due to these two states cancel each 
other and only the difference of the population proba- 
bility P(|fc| = 0.978) -P{]k\^ 0.982) contributes to the 
paramagnetic response. In Fig. 13 this corresponds to the 
difference between the filled mark and the open mark for 
each value of about ka/ir « 0.98. 

Figure 0fb) shows that the persistent current changes 
from the paramagnetic to the diamagnetic one as the 
strength of the impurities increases. When wq/^ = 2.0 in 
Fig. Etc), the difference P{\k\ = 0.978) - P{\k\ = 0.982) 
for n = AN becomes much smaller, whereas it is still 
very large for n = AN — 2 as shown in the inset of Fig. O 
Therefore, the persistent current /47V becomes weak com- 
pared to that of /4Ar-2 and, thus, the ensemble averaged 
persistent current shows the diamagnetic response as ob- 
served in the experimentsi^*^ A numerical simulation 
similar to the previous stud>«2a may be possible. Using 



the twisted boundary condition for the wave functions in 
an AB ring with a disordered zone inserted we think that 
it can confirm the present results. 

IV. SUMMARY 

In summary, we have presented a theoretical calcu- 
lation on the sign of ensemble averaged persistent cur- 
rent for the disordered AB rings. The experimentally 
observed diamagnetic response is attributed to the tran- 
sition process of the topmost electrons to the opposite 
branch for n = AN. It is shown that the backward scat- 
terings of the topmost electrons give rise to the diamag- 
netic currents in a natural way in agreement with the 
recent experimental observations. 
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